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ABSTRACT 

The accretion of satellites onto central galaxies along vast cosmic filaments is an 
apparent outcome of the anisotropic collapse of structure in our Universe. Numerical 
work (based on gravitational dynamics of iV-body simulations) indicates that satellites 
are beamed towards hosts along preferred directions imprinted by the velocity shear 
field. Here we use the Sloan Digital Sky Survey to observationally test this claim. We 
construct 3D filaments and sheets and examine the relative position of satellite galax¬ 
ies. A statistically significant alignment between satellite galaxy position and filament 
axis in observations is confirmed. We find a qualitatively compatible alignments by 
examining satellites and filaments similarly identified in the Millennium simulation, 
semi-analytical galaxy catalogue. We also examine the dependence of the alignment 
strength on galaxy properties such as colour, magnitude and (relative) satellite magni¬ 
tude, finding that the alignment is strongest for the reddest and brightest central and 
satellite galaxies. Our results confirm the theoretical picture and the role of the cosmic 
web in satellite accretion. Furthermore our results suggest that filaments identified on 
larger scales can be reflected in the positions of satellite galaxies that are quite close 
to their hosts. 

Key words: methods: data analysis - methods: statistical - galaxies: statistics - 
large-scale structure of Universe 


1 INTRODUCTION 

While the cosmic web (Joeveer, Einasto & Tago 1978; Bond, 
Kofman & Pogosyan 1996) is interesting in and of itself, its 
characteristics affect the haloes and galaxies that inhabit it. 
Observations and simulations have shown that the spatial 
distribution of satellite galaxies is not random, but rather is 
aligned with the major axes of host galaxies. Tormen (1997) 
was the first to suggest that galaxies fall into clusters/haloes 
anisotropically (see also Knebe et al. 2004; Aubert, Pichon 
& Colombi 2004). According to current knowledge the satel¬ 
lites are preferentially (not randomly) accreted along fila¬ 
ments and retain a memory of the large scales from which 
they came (e.g. Knebe et al. 2004; Libeskind et al. 2011). 

The alignment between the distribution of satellite 
galaxies and large-scale structure (described e.g. by the ve¬ 
locity shear tensor) is affected mainly by two dynamical 
processes: a possible pre-adjustment of satellites in the fila¬ 
ments, which points radially toward the host galaxy/group; 
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and the preferential accretion of satellites along those fila¬ 
ments. On smaller scales, the satellite distribution is also 
affected by the triaxial halo potential well. 

Numerous studies have characterized different types 
of alignments between host galaxies and their satellites. 
These include the alignment of the spatial distribution of 
satellites with the orientation of the host halo/galaxy (e.g. 
Knebe et al. 2004; Zentner et al. 2005; Kang et al. 2007; 
Azzaro et al. 2007; Libeskind et al. 2007; Wang et al. 
2008, 2010; Agustsson & Brainerd 2010; Wang et al. 2014b; 
Dong et al. 2014); the alignment of the orientation of satel¬ 
lites/subhaloes with respect to the centre of the host (e.g. 
Faltenbacher et al. 2008; Schneider et al. 2013); and the 
alignment between halo/group shape and the large-scale 
structure (e.g. Bailin et al. 2008; Paz et al. 2011; Zhang 
et al. 2013). Again, accretion along the filaments and the 
impact of tidal fields have been invoked to explain the ob¬ 
served alignment signals. All these studies reveal a signifi¬ 
cant alignment between galaxies and their surrounding envi¬ 
ronments. The alignment is observed to be stronger for red 
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galaxies/satellites and is almost absent between blue cen¬ 
trals and blue satellites. 

Despite the numerous studies that analyse the correla¬ 
tions between satellites and host galaxies/haloes, only re¬ 
cently, the alignment between satellites and surrounding 
large-scale structure have been analysed. Using A'^-body sim¬ 
ulations Libeskind et al. (2014) shows that the angular infall 
pattern of subhaloes is essentially driven by the shear tensor 
of the ambient velocity field. Dark matter (DM) subhaloes 
are found to be preferentially accreted along the direction 
of weakest collapse (es vector, the direction of filaments); a 
correlation is also found perpendicular to ei (sheet normal). 
Similar conclusions is reached in Forero-Romero & Gonzalez 
(2015) where the Local Group like pairs are analysed in sim¬ 
ulations. Also, the alignment of subhalo orbits is well corre¬ 
lated with the large-scale structure as defined by the velocity 
shear (Libeskind et al. 2012, 2013). 

In observations the distribution of satellites and veloc¬ 
ity shear are only analysed in nearby Universe. Lee, Rey & 
Kim (2014) show that the normal satellites and the dwarf 
satellites fall in to the Virgo cluster preferentially along the 
local filament and the local sheet, respectively. Libeskind 
et al. (2015) recently showed that planes of satellite galax¬ 
ies observed around the Milky Way (e.g. Lynden-Bell 1976, 
1982), the Andromeda galaxy (Ibata et al. 2013; Conn et al. 
2013) and Centaurus A (Tully et al. 2015) are well corre¬ 
lated with the velocity shear field, reconstructed by surveys 
of cosmic flows. Cautun et al. (2015) showed that the planes 
of satellite galaxies are generally in very good agreement 
with predictions from ACDM simulations, albeit this was 
challenged by Ibata et al. (2014). Pawlowski et al. (2012) ar¬ 
gued that filamentary accretion cannot explain the orbits of 
the Milky Way satellites. Also, Forero-Romero & Gonzalez 
(2015) conclude that the Andromeda-Milky Way system is 
in contradiction to the accretion along the filaments. Thus 
the alignment between substructures and their large-scale 
structure is of great interest in the study of the local Uni¬ 
verse. 

Analysing the spatial distribution of satellites in larger 
samples is generally more challenging because only few satel¬ 
lites are detected per primary galaxy, and because the three- 
dimensional location of the satellites is uncertain. This mea¬ 
surement requires that primaries and their satellites be 
stacked to obtain statistically viable samples. Owing to the 
advent of large galaxy surveys, a statistically robust treat¬ 
ment of satellite galaxies has become possible (e.g. Guo et al. 
2011, 2012; Wang & White 2012; Wang et al. 2014a). 

In the current work we use the Sloan Digital Sky Survey 
(SDSS) to analyse the alignment between satellite distribu¬ 
tion and the local filament/sheet orientation - indirect trac¬ 
ers of the local velocity shear tensor. For filament detection 
we use the Bisous model, which is based on a marked point 
process (Tempel et al. 2014c). We check our results against 
semi-analytical models of galaxy formation run on merger 
trees extracted from numerical simulations. Such alignment 
exist in V-body simulations (Forero-Romero & Gonzalez 
2015; Libeskind et al. 2014) but is never searched in a large 
observational datasets, and hence never checked directly for 
consistency with numerical work. 

The outline of the paper is following. In Sect. 2 we de¬ 
scribe the data and methods and in Sect. 3 we present our 
result. We conclude and discuss our results in Sect. 4. In Ap¬ 


pendix A using high-resolution V-body simulation we show 
that the Bisous model filaments and sheets are aligned with 
velocity shear tensor. 


2 DATA AND METHODS 

Throughout this paper we refer to isolated primary galax¬ 
ies (or just “primaries”) as the central galaxies that host 
systems of fainter satellites and fulfil a set isolation criteria 
as in Guo, Tempel & Libeskind (2015). A satellites posi¬ 
tion vector is simply the vector pointing from the primary 
to a given satellite. We also refer to the alignment between 
the anisotropic distribution of satellites with filaments axes 
simply as the alignment, if not specified otherwise. 

2.1 SDSS galaxies and filaments 

Galaxies and their satellites from both the spectroscopic and 
photometric samples in the SDSS data release 8 (York et al. 
2000; Aihara et al. 2011) are used in this study. The cat¬ 
alogue of isolated primaries and their satellites that are in 
filaments are the same as that in Guo et al. (2015). Pri¬ 
maries are selected from the spectroscopic sample, whereas 
satellite candidates can come from either the spectroscopic 
or photometric sample. We only consider isolated primaries 
in the sense that no other galaxy brighter than 0.5 mag less 
than the primary lies within a projected distance of 2Rinner 
(where Rinner is described below) and is sufficiently close in 
redshift. 

The limiting observed magnitude for satellite galaxies 
in the r filter is 20.5. For a satellite with a spectroscopic 
redshift to be assigned to a host, the difference in the red- 
shifts between the two must be less than 0.002. For satellites 
with only photometric redshifts, the difference between the 
photo -2 of the potential satellite and the primary’s spectro¬ 
scopic redshift of the host must be less than 2.5 times the 
photometric redshift error. All galaxies that are within a dis¬ 
tance Rinner/outer from a primary are considered satellites if 
they satisfy the redshift and magnitude criteria mentioned 
above. Because of the poorly constrained distances provided 
when only photometric redshifts are available, given a can¬ 
didate, it is not possible to determine if it is a genuine satel¬ 
lite or just a background/foreground galaxy (hereafter inter¬ 
loper) which is projected close to the primary. Statistically, 
genuine satellites tend to reside close to their primaries. As 
in Guo et al. (2015), we examine galaxies in the inner area 
(within a projected distance Rinner of a primary galaxy, here¬ 
after “inner” galaxies) and outer area (the annulus bounded 
by the projected distance Rinner and Router of the same pri¬ 
mary, hereafter “outer” galaxies). The values of Rinner and 
Router are the same as in Guo et al. (2015), which depend 
on the primary’s magnitude. We divided the primary sample 
into three r-band magnitude bins, each one magnitude wide 
and centred on Mr = —21.0, —22.0, and —23.0. The values of 
(Rinner, Router) are (0.3, 0.6), (0.4,0.8), and (0.55,0.9) Mpc, 
respectively. More details about primary and satellite galaxy 
selection are given in Guo et al. (2011, 2012). 

As shown in Guo et al. (2013, 2015), the number density 
of genuine satellites in the outer areas is usually lower than 
that in the inner area, since interloper galaxies dominate 
the outer areas. Also, interlopers are distributed randomly 
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with respect to the filaments. There should be no correlation 
between the angular distribution of interlopers and the lo¬ 
cal filaments axis. Therefore the alignment signals for outer 
galaxies must be diluted by these isotropically distributed 
interlopers. The details of measurement of the alignment 
signals and handling of interloper galaxies are described in 
Sect. 2.4. 

The catalogue of filaments is built by applying the 
Bisous process to the distribution of galaxies as described in 
Sect. 2.3. The method and parameters are exactly the same 
as in Tempel et al. (2014c). The assumed scale (radius) for 
the extracted filaments is 0.71 Mpc. Because the survey is 
flux-limited, the sample is very diluted farther away. Hence, 
we are only able to detect filaments in this scale up to the 
redshift 0.15 (« 640 Mpc). We refer to Tempel et al. (2014c) 
for the details of our filament finder. In addition to the fila¬ 
ment direction (es vector), we also estimate the orientation 
of sheets (ei and e 2 vectors) in the location of filaments as 
described in Sect. 2.3.1. 

In our analysis, we use only primary galaxies that are 
located in a filament, i.e. the primary distance from filament 
axis is less than 0.71 Mpc and the distance of the galaxy 
from the end point of filament (if the galaxy is outside a 
filamentary cylinder) is less than 0.14 Mpc. 

2.2 Semi-analytical galaxies and filaments 

The model galaxy sample was created using the publicly 
available Millennium Simulation (MS, Springel et al. 2005) 
to define the mass distribution and a semi-analytic model 
to place the galaxies within this density field. The simula¬ 
tion covers a large volume with a box length of 714 Mpc. 
The mass resolution used is ~ 10® Mq per particle. Galax¬ 
ies, whose properties are computed semi-analytically popu¬ 
late the DM haloes according the galaxy formation model 
GALFORM (Bower et al. 2006). It includes re-ionisation at 
high redshift, feedback from supernovae and stellar winds 
(in order to prevent the overproduction of the low luminos¬ 
ity galaxies), and outflows from active galactic nuclei (AGN) 
to suppress the formation of most luminous galaxies. The 
model was tuned to reproduce the observed A-band lumi¬ 
nosity at redshift 2 ; ~ 0. The low absolute magnitude limit 
of our model galaxy sample is about —16 in r band: this sets 
the magnitude limit for potential satellite galaxies. Fig. 5 of 
Guo et al. (2013) demonstrated that such resolution is good 
enough for a reliable comparison between the simulation’s 
satellite population and those in SDSS data. 

In the semi-analytical model, isolated primaries and 
their satellites are selected as described in sect. 3 of Guo 
et al. (2013), which is slightly different to the procedures 
applied to the SDSS light-cone survey data but with same 
physical motivation. The slight difference arises because the 
semi-analytical galaxies are in real-space while SDSS galax¬ 
ies are in redshift-space. For semi-analytical galaxies in mag¬ 
nitude bins, —21.0, —22.0, and —23.0, we define the term 
“inner” to mean within the sphere of the radius, 0.3, 0.4 and 
0.55 Mpc, respectively. The outer galaxies are those found 
in between ftnner and i?outer, where i7outer is 0.6, 0.8, and 
0.9 Mpc. The values of i7inner and 7?outer are the same as 
for SDSS galaxies, except they are measured in real space 
rather than in projection. The remaining isolation criteria 
for the semi-analytical galaxies are also the same as for SDSS 


galaxies. Since the inner and outer galaxies are selected in 
real space, there are no background/foreground interlopers 
among the satellite galaxies in this case. The simulation is 
thus used as a tool to identify the uncontaminated “true” 
signal. 

We have shown in Guo et al. (2013) that the estimated 
satellite luminosity functions and projected number density 
profiles from light-cone data (for SDSS sample, with back¬ 
ground subtraction) is not biased from those estimated in 
real space (for the model galaxy sample, without background 
subtraction). The estimated satellite luminosity functions 
and projected number densities are very similar. In this pa¬ 
per, we will directly compare the measured alignment sig¬ 
nals from light-cone data (for the SDSS sample) with the 
projected 3D position (for the semi-analytical galaxy sam¬ 
ple) as well. It is convenient to explore the dependence of 
the alignment signal on various properties for the simula¬ 
tion data, since there are (by constructions) no interlopers 
in the inner or outer areas. Because of the interlopers, the 
alignment signal for SDSS galaxies is necessarily weaker in 
most cases (e.g. see Agustsson & Brainerd 2010) compared 
with that from semi-analytical galaxies. That said, we do 
not explicitly attempt to recover the observational signal 
from the semi-analytical sample. Instead, we attempt to un¬ 
derstand the underlying dependencies of the alignment on 
galaxy properties. Therefore, in what follows we compare 
our results from the SDSS light-cone survey directly with 
those from semi-analytical galaxies, having in mind that we 
do not expect these two to be identical. 

The same filament finder is applied to both SDSS 
and semi-analytical galaxies to get the filaments catalogue. 
For the Bisous model, we use only galaxies brighter than 
— 19.0 mag in r band. In principle we could also use fainter 
galaxies to extract the filamentary network, however, we 
wish to avoid using satellites when extracting filaments. In 
this way, the alignment between satellite distribution and 
Bisous filaments and sheets is not built into the method 
(and hence trivial). 

For filament extraction, we use the same parameters 
for the Bisous model as we used for SDSS (see previous 
section). Primaries are classified as “in-a-filament” or “not- 
in-a-filament” in the same way for both samples. 

2.3 Bisous model for filament extraction 

We apply an object/marked point process with interactions 
(the Bisous process; Stoica, Gregori & Mateu 2005) to trace 
the filamentary network in the distribution of galaxies. This 
algorithm provides a quantitative classification which agrees 
with the visual impression of the cosmic web and is based on 
a robust and well-defined mathematical scheme. A thorough 
explanation of the method can be found in Stoica, Martinez 
& Saar (2007, 2010) and Tempel et al. (2014c); for conve¬ 
nience, a brief description is given below. 

The marked point process we use for filament detec¬ 
tion do not model galaxies/haloes, but models the structure 
outlined by galaxy/halo positions. This model approximates 
the filamentary network by a random configuration of small 
segments. We assume that locally galaxies may be grouped 
together inside a rather small cylinder, and such cylinders 
may combine to form a filament if neighbouring cylinders 
are aligned in similar directions. 
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The solution provided by our model is stochastic. 
Therefore, we find some variation in the detected patterns 
for different Markov chain Monte Carlo (MCMC) runs of 
the model. The main advantage of using such a stochastic 
approach is the ability to give simultaneous morphological 
and statistical characterization of the filamentary pattern. 

In practice, after fixing the approximate scale of the fil¬ 
aments, the algorithm returns the filament detection proba¬ 
bility field together with the filament orientation field. Based 
on these data, filament spines are extracted and a filament 
catalogue is built. Every filament in this catalogue is repre¬ 
sented as a spine; a set of points that define the axis of the 
filament. 

The spine detection we use is based on two ideas. First, 
filament spines are located at the highest density regions 
outlined by the filament probability maps. Second, in these 
regions of high probability for the filamentary network, the 
spines are oriented along the orientation field of the filamen¬ 
tary network. See Tempel et al. (2014c) for more details of 
the procedure. 

2.3.1 Complementing the Bisous filaments with sheet 
orientation 

As already mentioned, the Bisous model is probabilistic 
and gives the filament detection probability field (called 
the “visit map”). At the location of filament spines, the 
cross section of the visit map can be analysed. Assum¬ 
ing that filaments are located in sheets, the distribution of 
haloes/galaxies in filaments is elongated in the plane of the 
sheet. The real space elongation of haloes/galaxies will also 
affect the filament detection probability field. Hence, the 
cross section of visit map in the location of filaments should 
be elongated in the plane of the sheet. Consequently, the 
probabilistic nature of Bisous model can be exploited to de¬ 
fine the sheet normals in the location of filaments. This idea 
was already employed in Tempel & Libeskind (2013), where 
the orientation of galaxies with respect to the Bisous fila¬ 
ments and sheets were analysed. 

In Fig. 1 we show the examples of Bisous filaments de¬ 
fined using the visit map. In the right column the cross sec¬ 
tion of a filament is shown. The elongation of the visit map is 
modest, but this is expected since we use a cylindrical shapes 
to model the filamentary network and most importantly, ini¬ 
tially the Bisous model was designed to detect filaments not 
sheets. 

We designate the orientation of Bisous filaments as ea, 
whereas ea and ei are perpendicular to filaments and each 
other, and the latter defines the normal of the sheet defined 
using the Bisous model (see Fig. 1). 

As we show in Tempel et al. (2014b) and Appendix A, 
the Gi’s defined in the Bisous model are very well aligned 
with the eigenvectors of the velocity shear field, hence, the 
Bisous model can be used as an indirect tracer of the velocity 
shear. 

2.4 Measuring the alignment signal 

In order to examine (and quantify) if satellites are anisotrop- 
ically distributed with respect to filaments, the probability 
function P{0) = N{9)/{N{6)) (for both SDSS and semi- 
analytical galaxies) is measured. The alignment signal is 
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Figure 1. Illustration of Bisous model filaments and the large- 
scale structure axes defined by the Bisous model. Three exam¬ 
ples of filaments and their cross-sections are shown. Left panels 
show the filaments defined by the visit map. The filament direc¬ 
tion is indicated with es vector, whereas ei and e 2 vectors are 
perpendicular to it. Right panels show the cross-section of those 
filaments and the axes ei and e 2 . The axis ei is defined as the 
short semi-axis of the cross-section in the visit map. The axis e 2 is 
perpendicular to ei and es. The axis ei defines the sheet normal, 
where the filament is located. 

calculated using the kernel density estimation^. Here 6 is 
the angle between the filament axis and the satellite’s po¬ 
sition vector (relative to its host). For SDSS galaxies, only 
the component of 6 projected in the sky plane can be mea¬ 
sured. For the semi-analytical sample, we project the rela¬ 
tive satellite position onto an arbitrary plane (namely, one 
defined by two of the simulation’s box axes) and measure 
the angle 6 between the projected off-centre vectors and fil¬ 
ament axes. This is then compared with the measurements 
from the SDSS sample. N{9) and N{cos9) are the num¬ 
ber of the central-satellites pairs that subtend an angle 9 
from the local filamentary structure’s main defining axis. 
0 = 0° or cos 9 — 1 implies that the satellites are preferably 
distributed along the filament’s axis (hereafter, referred to 
as “alignment”), while 9 = 90° or cos 9 = 0 implies that 
satellites are preferably distributed perpendicular to this fil¬ 
amentary axes (hereafter termed “anti-alignment”). 

^ In Appendix A of Tempel et al. (2014a) we show that the kernel 
density estimation is better than a simple histogram and is more 
representative of the underlying probability distribution. 
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To calculate the alignment signal, we use all galaxies 
in the inner or outer regions. For semi-analytical galaxies, 
all galaxies in inner or outer regions are close to the central 
galaxy. For SDSS the inner and outer galaxy samples also 
include interloper galaxies that actually do not belong to 
the central galaxy. For the calculation of alignment signal, 
we do not try to remove interlopers or subtract any back¬ 
ground/foreground. Instead, we show the alignment signals 
for both inner and outer galaxies (using all galaxies in inner 
and outer regions)^. 

To estimate the significance of the alignment, we need 
to compare the measured probability function to the null- 
hypothesis of random (isotropic) angular satellite distribu¬ 
tion. However, due to selection effects, it is not granted that 
a random angular satellite distribution (nor the distribution 
of filaments axes with respect to the line-of-sight, Tempel, 
Stoica & Saar 2013) results in a uniform probability distribu¬ 
tion. Therefore we use a Monte Carlo (MC) approximation 
similar to the method used in Tempel et al. (2013) to esti¬ 
mate the case that the angular distribution of the satellites 
is random, regardless of filament axes and the confidence 
intervals for this estimate. This approach will take simul¬ 
taneously into account the biases in the filament detection 
and angular distribution of satellites. 

In order to do so, we generate 400 randomised samples 
(for the subsamples of relatively small volume, such as blue 
primaries, we generate 600 randomised samples) in which 
satellite position vectors with respect to their primaries are 
kept fixed, but the filament axes of the primary are shuf¬ 
fled and re-assign to each primary. The median and its 5% 
and 95% quantiles are then calculated by using these ran¬ 
dom samples. These confidence regions together with the 
Kolmogorov-Smirnov test (KS test) probabilities that the 
angle distributions are drawn from a uniform distribution 
are shown in the figures. 

3 RESULTS 

3.1 Satellite alignment in SDSS 

In Fig. 2, we show the probability distribution P{9) for the 
angle (in the sky plane perpendicular to the line of sight) 
between satellite position vectors and the filament/sheet 
axes (ei) for SDSS galaxies. Although the Rinner varies with 
the magnitude of primaries, for the angle 9 we can mea¬ 
sure the mean excess probability distribution P{9) for all 
primaries (including all galaxies in three magnitude bins 
Mp = —21, —22, —23) by simply stacking them together. 
We will explore the possible dependence of the excess prob¬ 
ability on the properties of galaxies later. In each panel of 
Fig. 2, we show the average angle (0i/o) and the Kolmogorov- 
Smirnov (KS) test probability that the sample is drawn 
from the randomised distribution for inner/outer galaxies. 

Fig. 2 shows that the satellites around primaries are 
preferentially distributed along the filament axis (ea vector) 

^ Note that interlopers mostly affect the strength of the align¬ 
ment signal. In general, if the strength of the alignment signal 
is important, the background/foreground should be statistically 
taken into account (e.g. see Wang & White 2012; Cautun et al. 
2015). 


and have anti-alignment with the sheet normal (ei vector) 
where the filament is located. This results is in accordance 
with the predictions from pure A^-body simulation (Libe- 
skind et al. 2014). Similarly in simulations, the alignment 
between satellite distribution and ei, es vector is statisti¬ 
cally significant (see the KS test values), while the alignment 
between satellites and 02 vector is very weak or absent. 

In order to compare our results with simulations, in 
the next section we calculate the same alignments for the 
galaxies in the Millennium simulation. 


3.2 Satellite alignment in Millennium simulation 

In this section we examine satellite positions with respect to 
the ei vectors as defined by Bisous model using the Millen¬ 
nium simulation. Fig. 3 shows the alignment between satel¬ 
lites and ei vectors in real space (upper row) and in pro¬ 
jection (lower row). In both cases a statistically significant 
alignment is seen, and in qualitative agreement with the 
alignment measured for the observed SDSS galaxies. The 
strength of the alignment between satellite positions and 
the ei vectors is (obviously) weakened when computed in 
projection; it is however not erased. 

For both SDSS and semi-analytical primaries, the align¬ 
ment results show similar trends that the satellites around 
them preferentially distribute along the filaments axes and 
have anti-alignment with the sheet normal. For SDSS galax¬ 
ies, this alignment is stronger for inner galaxies than for 
outer galaxies (see Fig. 2), which is opposite for semi- 
analytical galaxies (Fig. 3). The signals of alignment for 
outer galaxies in SDSS sample could be diluted because most 
galaxies in the outer areas are interlopers whose distribution 
is isotropic with respect to the filament axis. However, for 
the semi-analytical galaxy sample, the outer galaxies are ac¬ 
tually relevant to the primaries - they are not far from the 
primaries in distance. Furthermore, the distribution of these 
outer galaxies may be influenced directly by the filaments. It 
is possible that these outer galaxies are themselves also pref¬ 
erentially distributed along the hlament axis and are even 
more strongly aligned with it. This is because the distribu¬ 
tion of outer galaxies is less affected by the primary’s DM 
halo (compared with inner galaxies that may well be orbiting 
within the host’s DM halo). After galaxies are accreted onto 
the main DM halo, the “memory” of the alignment with the 
Hlament axis may be forgotten or weakened (e.g. see Vera- 
Ciro et al. 2011). This interpretation is in accordance with 
the alignment of galaxy pairs in filaments. Namely, Tempel 
& Tamm (2015) shows that the alignment of loose galaxy 
pairs with filaments is much stronger than that for close 
galaxy pairs. 

Recall that the strength of the alignment measured from 
the SDSS and Millennium simulation sample is expected 
to be different (see Sect. 2.4). The strength of alignment 
depends on several factors, (examined below) including the 
number of background/foreground galaxies in the inner or 
outer areas, the properties of primaries and satellites used 
for the measurement. 
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Figure 2. The alignment of satellites in SDSS as a function of the angle 6 (measured in the plane perpendicular to the line of sight) 
between the satellite position vector (relative to the primary galaxy) and the orientation of structures (e;) as defined in the Bisous model. 
The vector 03 (right panel) designates the filament orientation, while ei (left panel) gives the normal of the sheet where the filament is 
located. The alignment is shown for inner (close to the primary) and outer (farther away from the primary) satellites. The definition for 
inner and outer satellites is given in the text. The filled regions show the 5th and 95th percentile spread for a randomised distribution. 
On each panel, we also give the average angle and the Kolmogorov-Smirnov (KS) test probability that the angle distribution is drawn 
from a randomised distribution. Satellite galaxies in the SDSS are statistically well aligned with the filamentary structure 03 . 


:' ' ' ' 1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' '; 
r Simulation (real): ei -■ 

1.15 

1.10 

:' ' ' ' 1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' ': 
r Simulation (real): e 2 ^ 

1.15 

1.10 

:' ' ' ' 1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' '; 
r Simulation (real): 63 /-z 

^ yr 

- V ^ 

1.00 

0.95 

: N : 

1.00 

0.95 


r 


- - 

r • — Outer: <008 00 > = 0.48 Pks < 

ri ■ ■ ■ 1 ■ ■ ■ ■ 1 ■ ■ ■ ■ 1 ■ ■ ■ ■ 1 ■ ■ ■ 

0.90 

0.85 

r • — Outer: <cos 0o > = 0.50 p^g = 7.7 -1 0 ^ 

ri ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ 

0.90 

0.85 

i- • — Outer: <cos0o> = O.52 p^g<1.0 10 ® 

i-i ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ 



0 20 40 60 80 0 20 40 60 80 0 20 40 60 80 

see 


Figure 3. The real (upper row) and projected (lower row) space alignment between anisotropic distribution of satellites and e; vectors 
in Millennium simulation. The definition of is given in Sect. 2.3.1; 03 gives the filament orientation and ei defines the normal of the 
sheet where the filament is located. The designation of lines and labels are the same as in Fig. 2. 


3.3 Dependence on galaxy properties 

In this section, we split our primaries into different sub¬ 
samples. We then explore the possible dependence of the 
alignments on the properties of primaries or satellites. 


3.3.1 Colour dependence of alignment signal 

Firstly, we split the primaries into subsamples by ° °((/ — r) 
colour of the primaries. We use two different lines, ^'^(g — 
r) = 0.15 - 0.024 Mr and ° °{g - r) = -0.28 - 0.04 Mr 
to divide the primaries into red and blue subsamples for 
SDSS and semi-analytical galaxies, respectively, due to the 
distribution of SDSS galaxies in the colour-magnitude di¬ 
agram is shifted along the g — r axis compared to that of 
the semi-analytical galaxies (refer to Guo et al. 2013, for 
the choice of the lines). Fig. 4 shows the measured align¬ 
ment for red and blue primaries respectively. For the red 


primaries in both samples, the strength of alignment of inner 
and outer galaxies are both very pronounced. The KS test 
probability for distribution of inner galaxies around SDSS 
red primaries is as low as ~ 10“^. However, for SDSS blue 
primaries, the result of angular distribution of inner galaxies 
and outer galaxies with respect to filament axes is statisti¬ 
cally isotropic (Fig. 4, bottom left panel), which is different 
to the results from the simulation (Fig. 4, panel on the bot¬ 
tom row, third column). For model blue primaries, the inner 
galaxies around them show a weak preference of distribut¬ 
ing along the filament axes, while a strong preference of 
distributing along the filament axes for outer galaxies. This 
suggest that the alignment for outer galaxies might be weak¬ 
ened when they fall into the main halo (and became inner 
galaxies). In SDSS the alignment is not observed, probably 
because of interloper galaxies in the sky projection. 

In Fig. 4 we show the alignment for all vectors and 
the dependence on colour is the same for all of them. Since 
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Figure 4. The alignment between anisotropic distribution of satellites and vectors in SDSS (two left-hand columns) and in Millennium 
simulation (two right-hand columns). The alignment is shown for red and blue primaries. The designation of lines and labels are the 
same as in Fig. 2. 
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Figure 5. The alignment between anisotropic distribution of 
satellites and filament orientation (03 vector) in SDSS (upper 
row) and in Millennium simulation (lower row). The alignment is 
shown for bright (Mp = —23) and faint (Mp = —21) primaries. 
The designation of lines and labels are the same as in Fig. 2. 

the dependence is strongest for es vector, we continue by 
presenting only this alignment. The alignment for ei is sim¬ 
ilar, but weaker and omitted in the interest of clarity. 

3.3.2 Dependence on primary magnitude 

We split the primaries into subsamples by magnitude and 
explore the dependence of the alignment on this quantity. 
Fig. 5 shows the resulting alignments for these subsamples. 
For both SDSS and model galaxy subsamples, the strength 
of alignment slightly depends on the luminosity of primaries; 
although for SDSS galaxies, different number of interlopers 
could also cause the variation in the strength. The (0i) for 


Figure 6. The alignment between anisotropic distribution of 
satellites and filament orientation (03 vector) in SDSS (upper 
row) and in Millennium simulation (lower row). The alignment is 
shown for bright and faint (relative to the luminosity of primary) 
satellites. The designation of lines and labels are the same as in 
Fig. 2. 


inner galaxies decrease from 44.70 to 44.42 from primaries 
of low luminosity to high luminosity. The {9{) for primaries 
in magnitude bin Mp = —22 is 44.70, which is roughly the 
same as that for magnitude bin Mp = —21. However the 
change of overall shape of distribution P{6) with the pri¬ 
mary luminosity is consistent with our previous statement. 
Moreover, the results from semi-analytical galaxies, which 
is not affected by interlopers, also show similar dependence: 
the strength of the alignment will increase with the luminos¬ 
ity of primaries. The luminosities of primaries are loosely 
correlated with their masses. Therefore, we could also say 
that the strength of alignment depends on the mass of the 
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primary galaxy (analogous to the trend found for DM halo 
shapes in Libeskind et al. 2013). 

3.3.3 Dependence on satellite magnitude 

Here, we split the satellites around the all primaries into 
subsamples by the luminosity of satellites, bright satellites 
Afsat < Mp + 2, middle luminous satellites Mp + 2 < Msat < 
Mp + 3, and faint satellites Mp + 3 < Msat < Mp + 4. 
Then we explore how the alignment depends on the rela¬ 
tive satellite-host luminosity. Fig. 6 shows the alignment for 
satellites in bright and faint subsamples. Again, both results 
from SDSS and model samples imply that the strength of 
alignment slightly depends on the luminosity of satellites. 
The {6i) for three SDSS subsamples (model subsamples) are 
44.06, 44.38 and 44.70 (42.37, 42.55, 42.91), which is anti¬ 
correlated with satellites luminosity. The trend for {6o} is 
the same in both cases. This shows that the alignments for 
bright satellites are stronger. The results also indicate that 
the measured alignment with respect to is not affected 
by the faint magnitude limit in the SDSS and in simula¬ 
tions: the alignment is qualitatively the same regardless of 
the luminosity of satellites. 

3.4 Intrinsics of the alignment signal 

At this point in our study it is relevant to highlight that 
numerous observational studies (since the seminal work of 
Holmberg 1969) have examined the angular distribution of 
satellites around primaries (Zaritsky et al. 1997; Brainerd 
2005; Agustsson & Brainerd 2006, although not all of these 
studies are in full agreement). Perhaps the most relevant 
work in this context is Yang et al. (2006) who found that 
satellites in the SDSS tend to be found along the main axis 
of their host galaxies. We find similar alignment (see Ap¬ 
pendix B2). On the other hand, the major axis of galaxies 
was also found to be correlated with the axes of filaments 
(e.g. Tempel et al. 2013), which is again corroborated here 
in Appendix Bl. Therefore we wish to explore if the align¬ 
ment between the anisotropic distribution of satellites and 
hlament axes shown above found here is intrinsic or just the 
result of the collective effect of two former alignments that 
are already well known. Hereafter, the known alignments are 
referred as “established” alignments. 

In order to check this we collect all the angles between 
the filament axes and primary galaxy major axes in the plane 
of the sky. Then for each primary, we randomly select a new 
angle from this distribution. We thus generate an artificial 
filament axes for each primary according to this new picked 
angle^. The resulting alignment signal between these artifi¬ 
cial filament axes and the primary major axes is by construc¬ 
tion statistically identical to that with the real filaments as 
shown in the top panel of Fig. 7. A crucial aspect of this 

® An alternative possibility is to randomise the primary galaxy 
major axes with respect to filaments, while fixing the filament 
orientation and satellite distribution. In this case, the alignment 
between artificial major axes and satellite positions should be 
analysed. Statistically speaking, these two approaches are iden¬ 
tical in the ideal case that all satellite are aligned with filament 
axes (since in both cases same angles are randomised). 



0 10 20 30 40 50 60 70 80 90 

0 


Figure 7. Test of the intrinsics of the alignment between the fil¬ 
ament axes and anisotropic distribution of satellites. Top panel: 
the alignment between real/artificial filament axes and the ma¬ 
jor axis of red primaries. Bottom panel: the alignments between 
anisotropic distribution of satellites with artificial filament axes 
(thick red and purple lines) and with real filament axes (thin blue 
and green lines). 


test is that it leaves the alignment between the satellite po¬ 
sition and galaxy major axis, undistributed. Therefore any 
direct correlation between the filaments and the anisotropic 
distribution of satellites is randomised, and the remaining 
correlation is due to the combined effect of two established 
alignments. If the previous alignments between real filament 
axes and satellite positions is intrinsic, the alignment signal 
will be much weaker when artificial filaments axes are con¬ 
sidered, since the direct correlation is broken. If the satellite 
position - filament axis alignment is just the consequence of 
the two established alignments mentioned earlier, a similar 
alignment signal should also exist when the artificial fila¬ 
ments are considered, since these two established alignments 
are preserved and statistically same as for real filaments. 

Suppose satellite positions are locked to the major axis 
of their host galaxy. In this case randomising the angle be¬ 
tween the filament and the host major axis will have the ef¬ 
fect of randomising the angle between the filament and the 
satellite distribution. In this case however, the distribution 
of angles between the satellite galaxy position and the fila¬ 
ment will not change since by construction the distribution 
of angles between the galaxy and filament has not changed. 
Now if on the other hand, satellite galaxies are locked to 
the filaments, then randomising the angle between filament 
and host major axis will have (1) no effect on the satel¬ 
lite galaxy position-host major axis alignment (this has not 
been touched) and (2) no effect on the galaxy major axis- 
filament alignment (by construction). The effect will only be 
to erase the direct alignment between the satellites and the 
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Figure 8. The alignment between anisotropic distribution of 
satellites and filament orientation (03 vector) in SDSS (upper 
row) and in Millennium simulation (lower row). The alignment is 
shown for primary galaxies that are in inner part of filaments or 
in outer part of filaments. The designation of lines and labels are 
the same as in Fig. 2. 

filament. The remaining alignment is the collective effect of 
two established alignments. 

The result of alignment between the artificial filament 
axes and the major axis of red primaries is shown in the 
bottom panel of Fig. 7. It shows there is statistically no 
alignment (or the alignment is very weak) between the arti¬ 
ficial filament axes and anisotropic distribution of satellites, 
which implies that the well established alignments are not 
sufficient to explain the measured alignment between the 
anisotropic distribution of satellites and filament axes. In 
other words: the alignment of satellites with the filamentary 
structure is intrinsic. 

As one final check, we split our primaries according to 
their distance to the filament axis and examine the satellite 
alignment, since Tempel & Libeskind (2013) also found the 
dependence of this property on major axis-filament align¬ 
ment: the alignment is stronger in outer parts of filaments. 
The results are shown in Fig. 8 . For galaxies further away 
from the filament axis, the satellite-filament alignment is 
weaker (even for red primaries). This is another indicator 
that the alignment we found is not trivial results of combined 
effects, since otherwise the alignment should be stronger in 
outer parts of filaments (since the alignment between galaxy 
major axis and filaments is stronger in outer parts). 


4 DISCUSSION AND CONCLUSIONS 

In the present paper, we extended the probabilistic model for 
filament detection (called Bisous model; Tempel et al. 2014c) 
and showed that the stochastic nature of Bisous model can 
be used to trace the orientation of large-scale structure in 
the location of galactic filaments. We analysed the alignment 
between satellite distribution around host galaxies and cos¬ 
mic filaments in the Millennium simulation and in SDSS 
observations. In observations, results can only be obtained 
in redshift space; these are (directly) compared with the 
simulation results, which are computed in real space. When 
examining the semi-analytical galaxies we have the benefit 


of omitting interlopers, thus making any comparison more 
qualitative than quantitative. We expect that interlopers 
and nuances related to the association of satellites with pri¬ 
maries in the observations will have the effect of weakening 
any inherent alignment. This will not affect our ability to 
quantify the dependence of the alignments on various prop¬ 
erties of the primary or satellite galaxies. 

The probabilistic nature of Bisous model allows us to 
define the sheets where the galactic filaments are located. 
The sheet normal, designated as ei, is well aligned with 
the minor principal axis of the local velocity shear (See Ap¬ 
pendix A). This extends the work by Tempel et al. (2014b), 
where the orientation of Bisous model filaments were com¬ 
pared with the major principal axis ( 03 ) of the velocity 
shear. It implies that in observations the Bisous model can 
be used as an indirect tracer of velocity shear tensor in lo¬ 
cation of galactic filaments. Our main results can be sum¬ 
marised as following. 

• Using the SDSS observations, we found weak but sta¬ 
tistically significant signal of alignment between the angular 
distribution of satellite galaxies around the isolated primary 
galaxies in filaments and the direction of filaments/sheets 
where those primaries are located. These findings confirm 
the predictions from simulations (Libeskind et al. 2014; 
Forero-Romero & Gonzalez 2015). 

• Using the Millennium simulation populated with semi- 
analytical galaxies, we showed that both observations and 
simulations show a qualitatively compatible alignment be¬ 
tween satellite galaxies and the orientation of filaments. 
These results are consistent and complement the predictions 
from pure A-body simulation (Libeskind et al. 2014). 

• The alignment signal is much stronger for red primary 
galaxies and is very weak or absent for blue primaries. This 
trend is visible both in observations and in simulations. We 
also see a weak trend that the alignment signal is stronger 
for brighter primaries and for brighter satellite galaxies. 

Our results confirm the claim that the velocity shear 
field dictated the infall of satellites into host systems (Libe¬ 
skind et al. 2014; Forero-Romero & Gonzalez 2015). A-body 
simulations show that satellites are aligned with the princi¬ 
pal axis corresponding to slowest collapse ( 03 ) of the ve¬ 
locity shear and show an anti-alignment with the principal 
axis corresponding to fastest collapse (ei), whereas there 
are almost no alignments with 62 vector (intermediate prin¬ 
cipal axis). This work confirm those DM only studies using 
the Millennium simulation semi-analytical galaxies as well 
as SDSS observations. 

Libeskind et al. (2014) show that the infall of satellites 
is universal in the sense that the strength of satellite beam¬ 
ing does not depend on the redshift or mass of the halo 
or satellite. In the present study we show that the satel¬ 
lite alignment depends on the colour/luminosity of the pri¬ 
mary galaxy: brighter satellites and brighter hosts display 
a stronger alignment. Indeed this is in perfect agreement 
with the numerical studies of Libeskind et al. (2014), which 
showed a similar mass dependence of the beaming effect. 
However, for blue/faint primaries, the satellite distribution 
around them may not reflect the infall direction of satellites, 
or the related information are disrupted when satellites fall 
into the main system. Agustsson & Brainerd (2010) found 
the satellites distribution with respect to the major axis 
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Figure 9. The alignment between anisotropic distribution of 
satellites and filament orientation in SDSS (upper row) and 
in Millennium simulation (lower row) for a fixed stellar mass 
(10 < log(Mgteiiar) < H)- The alignment is shown for red (right 
column) and blue (left column) primary galaxies. The designation 
of lines and labels are the same as in Fig. 2. 


of primary galaxies have similar dependence on the colour 
of primary galaxies. However they found that the isotropic 
satellite distribution around blue primaries are due to that 
the satellite distribution dependent also on the projected 
distance of satellites to the primaries. The satellites close 
the primaries actually preferably distribute along the major 
axes of primaries, while those far from primaries distribute 
randomly or preferably along the minor axes. We also tried 
to explore the possible dependence of anisotropic distribu¬ 
tion on the projected distance to the primaries for using both 
SDSS and model galaxies. However we found no statistical 
evidence supporting that. To give a conclusive answer to 
this question, a more detailed study based on simulations is 
required. 

The alignment dependence on primary colour is an in¬ 
teresting result and requires more attention. Logically, the 
colour dependence may simply reflect the mass (or lumi¬ 
nosity) dependence of the alignment: redder primaries are 
also more massive and more massive galaxies show stronger 
alignment signal. To investigate this effect, we selected pri¬ 
maries on a fixed stellar mass^ bin (10 < log(Mateiiar) < 11 ) 
and calculated the alignment signal for blue and red pri¬ 
maries®. The resulting alignment signals for SDSS and sim¬ 
ulations are shown in Fig. 9, where we see that the depen¬ 
dence on primary colour is intrinsic. However, the number of 
satellites (and significance) for blue primaries is lower than 
for red primaries, because red primaries have more satellites 
than blue ones of the same stellar mass (Wang & White 
2012 ). 

Since redder colour typically indicate older stellar pop¬ 
ulation, these results suggest that a significant alignment 


^ For SDSS stellar masses were estimated using the KCORRECT 
algorithm (Blanton &; Roweis 2007). For semi-analytical galaxies, 
the stellar masses are direct outcome from GALFORM. 

® We also conducted the analysis using narrower stellar mass 
bins. The conclusions in this case remain the same, however, the 
significance is reduced since the number of primaries decreases. 


only exists in haloes with a relatively old stellar popu¬ 
lation. Agustsson & Brainerd (2010) argue that satellites 
that recently enter to the host halo show considerably less 
anisotropy than do those which entered their host’s halo at 
earlier times. Perhaps the trajectories of satellites are af¬ 
fected by the potential well of the surrounding filamentary 
environment. 

Clearly, the satellite-filament correlation must hold 
some interesting clues regarding galaxy formation. The mea¬ 
sured alignments between satellite positions and the sur¬ 
rounding filamentary environment may indicate that galaxy 
formation is affected by large-scale environment. The align¬ 
ment signal may be natural result of how satellites accrete 
via streams along the direction of the filaments. The results 
also suggest that the red central galaxies are significantly 
affected by such accretion. 

Recently, Lee & Choi (2015) carried out similar study 
as presented in this work. Lee & Choi (2015) used the SDSS 
DR7 observations and reconstructed the velocity shear field 
using linear approximation. Motivated by the work of Libe- 
skind et al. (2014), they show that satellites around isolated 
galaxies are aligned with 63 and show anti-alignment with 
ei. They conclude that the alignment with respect to veloc¬ 
ity shear is universal as predicted in Libeskind et al. (2014). 
In general, our study is in good agreement with the work 
by Lee & Choi (2015). There are two important differences. 
First of all we see a dependence on primary galaxy colour 
and luminosity. Secondly we use a fully independent method 
(the Bisous model) to compute filamentary structure. The 
agreement with Lee & Choi (2015) supports the notion that 
satellites are not only beamed towards their hosts by the 
cosmic web, but that this effect can be seen long after the 
accretion process has ended. 

As a complementary study, Tempel & Tamm (2015) 
shows that galaxy pairs in SDSS are also very strongly 
aligned with galactic filaments where they are located. 
Adding the results from this paper, we can conclude that 
not only satellite galaxies, but also galaxy pairs and po¬ 
tentially all structures are elongated along the galactic fila¬ 
ments. These results will put important constraints on the 
formation and evolution of galaxies and groups in the cosmic 
web. 
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Figure Al. The alignment between Bisous model filaments and sheets (e^ Bisous vectors) and the underlying velocity field as defined by 
the principal axes of velocity shear tensor (e; vectors). Left panel shows the alignment between ei Bisous and e^; middle panel between 
62. Bisous and 6;; right panel between 63 Bisous and e^. The vectors 63 Bisous and 63 define the orientation of filaments, while 63 Bisous 
and 61 define the sheet normals. With dotted grey lines we show the random distribution between 61.2,Bisous and 62,1 assuming perfect 
alignment between 63. Bisous and 63. Strong alignment between the same principal axes of Bisous structures and shear tensor is visible, 
hence, the Bisous model structures trace the underlying velocity field. In lower row, the orientations of 61.2,Bisous are randomised, while 
keeping the 63 Bisous orientation. Comparing with upper row we note that the real orientation of 63 Bisous and 62 Bisous is not drawn 
from a random distribution. 


APPENDIX A: CORRELATION BETWEEN 
VELOCITY SHEAR TENSOR AND BISOUS 
MODEL FILAMENTS AND SHEETS 

Here we analyse the correlation between velocity field (shear 
tensor) and the structures defined in the Bisous model using 
a pure DM N-body simulation. We extend the work pre¬ 
sented in Tempel et al. (2014b), that showed that the direc¬ 
tion of filaments defined using point processes applied to the 
halo distribution in DM N-body simulation matches very 
well the eigenvector of the velocity shear tensor correspond¬ 
ing to slowest collapse. In Tempel et al. (2014b) only the 
orientation of filaments were compared with the 63 vector 
of the local velocity shear. Now we employ the probabilistic 
nature of the Bisous process and define the galactic sheets 
(their normals) in the location of filaments (see Sect. 2.3.1). 
The aim is to test whether the halo/galaxy distribution in 
filaments can be used as a tracer of velocity shear tensor. 
For that the sheet normal defined in the Bisous model are 
compared with the eigenvectors of the local velocity shear 
tensor. 


Al V-body simulation and the velocity shear 
tensor 

The simulation we use in this Appendix is exactly the same 
as used in Tempel et al. (2014b). For convenience a brief 
description is given below. 

A DM-only Y-body simulation is run assuming the 
standard ACDM concordance cosmology, in particular a 
flat universe with cosmological constant density parameter 
Qa = 0.72, matter density parameter flm = 0.28, a Hub¬ 
ble constant parameterized by Ho = 100 h km s“^Mpc“^ 
(with h = 0.7), a spectral index of primordial density fluc¬ 


tuations given by ris = 0.96, and mass fluctuations given by 
as = 0.817. 

The simulations span a box of side length 64 /i“^Mpc 
with 1024^ particles, achieving a mass resolution of ~ 1.89 x 
10^ h~^MQ and a spatial resolution of 1 ft“^kpc. The pub¬ 
licly available halo finder AHF (Knollmann & Knebe 2009) 
is run on the particle distribution to obtain a halo catalogue. 
Only haloes more massive than 10® h~^MQ are considered 
in this work. 

We quantify the cosmic web by means of the velocity 
shear tensor. This method is described in detail in Hoffman 
et al. (2012). The salient aspects are highlighted here, in 
brief. The cosmic velocity field is calculated using a “Clouds- 
in-Cell” (CIC) algorithm on a 256® grid. The velocity (and 
density) fields are then smoothed with a Gaussian kernel 
equal to at least one cell (i.e. ramooth = 0.25 fi“®Mpc) in 
order to get rid of the spurious artificial cartesian grid intro¬ 
duced by the CIC. In practice the smoothing sets the scale 
of the calculation. In current study, we use only the smooth¬ 
ing scale 0.5 /i“®Mpc. The velocity shear tensor is defined as 

~ calculated by means of fast 

Fourier transform (FFT) in fc-space. The velocity shear ten¬ 
sor is then diagonalized and its eigenvectors and eigenvalues 
are identified. We denote the eigenvectors as 63 , 62 and 63 , 
where 63 is the direction of fastest collapse and 63 gives the 
direction of slowest collapse (e.g. the direction of filaments). 

Note that the velocity shear field is identical to the tidal 
field, defined as the Hessian of the potential, namely Tap = 

Q 2 j 

dr drp ’ "'hen smoothed on large enough (i.e. > few Mpc) 
scales. 
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A2 Orientation of Bisous filaments and sheets 
with respect to the underlying velocity field 

Tempel et al. (2014b) showed that the Bisous filaments are 
very well aligned with the slowest collapse in underlying ve¬ 
locity field. Here, we extend the study and analyse the ori¬ 
entation of Bisous model filaments and sheets with respect 
to the all three eigenvectors of the velocity shear tensor. 
For that we calculate the distribution of cosine of the an¬ 
gles between the Ci,Bisous (Bisous model) and Ci (velocity 
shear) vectors. The angle between Bisous structures and ve¬ 
locity field is found in the location of filament axes. For every 
point in the filament spine, the angle is found between the 
orientation of that filament point (ci,Bisous) and the veloc¬ 
ity shear eigenvector (e^) defined in the closest cell. Since 
the scale and spacing of points along the filament spine is 
0.5 /i“^Mpc, which is larger than the cell size (0.25 /i~^Mpc), 
interpolation of velocity shear to the filament axis is not re¬ 
quired. 

Fig. Al shows the alignment between the structure ori¬ 
entations in the Bisous model (ei, Bisous, 62,Bisous, 63,Bisous) 
and the velocity shear (ei, 62 , ea). Fig. Al shows that in 
addition to the very strong alignment between ea,Bisous and 
ea, there is also a strong alignment between ei,Bisous and 
Cl, and between 62 ,Bisous and ea. This shows that the sheet 
normal ei,Bisous defined in the location of Bisous filaments 
is aligned with the fastest collapse in the underlying veloc¬ 
ity field, e.g. the Bisous model can be used as an indirect 
tracer of underlying velocity field in cases where the velocity 
information is not available (e.g. in observations). 

It is expected that if ea^isous and ea are perfectly 
aligned, then even for completely random distribution for 
ei, Bisous, the latter is somewhat aligned with ei and 02 
in real space. To show that this is not so in Bisous struc¬ 
tures, we randomised the ei^Bisous while keeping the ea^isous 
orientation (direction of filament). The alignments for ran¬ 
domised ei, 2 ,Bisous are shown in lower panel of Fig. Al. If 
we compare the alignments in upper and lower panels in 
Fig. Al, we clearly see that the Bisous model ei, 2 ,Bisous vec¬ 
tors are not drawn from random distribution. The alignment 
between the sheet normal in Bisous model (ei,Bisous) and ve¬ 
locity field eigenvector ei is real. 


APPENDIX B: KNOWN ALIGNMENTS OF 
SATELLITES AND PRIMARIES 

B1 The alignment between the major axis of 
primaries and filament axes 

In this section, we corroborate that the major axis of pri¬ 
maries aligns with the filament axes, which has already been 
studied (e.g. Tempel et al. 2013). Fig. B1 shows the excess 
probability of primaries as a function of the angle between 
the major axis of primaries and filaments axes. The results 
corroborate that for red primaries, the major axis preferen¬ 
tially align with the filaments axes. Here we used the de Vau- 
couleurs fit position angle {<j)) as the proxy of the major axis 
of the primaries. 
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Figure Bl. The excess probability of primaries as a function 
of angles between the major axis of primaries and the filament 
axes. The results for red and blue primaries are shown in left and 
right panel, respectively. Dotted lines with filled region show the 
null-hypothesis together with its 95% confidence limit. 



Figure B2. The normalised number of satellites as a function of 
the angle between the satellite-primary position vector and the 
main axis of the primary galaxies (position angle cj>) in sky plane. 
The alignment is shown for blue (top row) and red (bottom row) 
primary galaxies, and for primaries in filaments (left column) and 
not in filaments (right column). 


B2 The alignment between anisotropic 

distribution of satellites and the major axis of 
primaries 

In this section, we corroborate that anisotropic distribution 
of satellites align with the major axis of primaries both 
in-filaments and not-in filaments. The result are shown in 
Fig. B2, which indicate that for primaries both in-filaments 
and not-in-filaments, the anisotropic distribution of satel¬ 
lites align with the major axis of red ones, but randomly 
distributed with respect to the the major axes of blue ones. 
It is consistent with the results from other studies (e.g. Yang 
et al. 2006). 
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